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Critical properties of spherically symmetric accretion in a fractal 
medium 
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ABSTRACT 

Spherically symmetric transonic accretion of a fractal medium has been studied in both the 
stationary and the dynamic regimes. The stationary transonic solution is greatly sensitive to 
infinitesimal deviations in the outer boundary condition, but the flow becomes transonic and 
stable, when its evolution is followed through time. The evolution towards transonicity is more 
pronounced for a fractal medium than what is it for a continuum. The dynamic approach also 
shows that there is a remarkable closeness between an equation of motion for a perturbation 
in the flow, and the metric of an analogue acoustic black hole. The stationary inflow solutions 
of a fractal medium are as much stable under the influence of linearised perturbations, as they 
are for the fluid continuum. 
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1 INTRODUCTION 

Accretion processes involve the non-self-gravitating flow dynamics 
of astrophysical matter under the external gravitational influence of 
a mass ive astrophysica l object, like an ordinary star or a compact 
object dFrank et alj|2002h . A paradigmatic model of an astrophysi- 
cal accreting system is that of spherically symmetric infall on to a 
centra l accretor. Ever since the seminal work published by iBondJ 
(1952), which effectively launched the subject in the form in which 
it is recognised today, the problem of spherically symmetri c flows 
has been revisited time and again from various angles 
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Roy 12007 ). This abiding appeal of the spherically symmetric 



model is explained by the fact that almost always it lends itself to 
an exact mathematical analysis, and in the process it allows a very 
clear insight to be had into the underlying physical principles. 

Ease of mathematical manipulations, however, is not the only 
reason why spherically symmetric flows are regularly invoked in 
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accretion-related literature. The details of the physics of many as- 
trophysical flows are, in fact, very faithfully described and under- 
stood with the help of this relatively simple model. Accretion of the 
interstellar medium (ISM) is a case in point. 

While formal fluid dynamical equations in the Newtonian con- 
struct of space and time — which would involve a momentum bal- 
ance equation (with gravity as an external driving force), the con- 
tinuity equation and an equation of state — suffice to a great ex- 
tent in shedding light on the accretion of the ISM, it must at the 
same time be recognised that the ISM is not entirely to be seen as 
a fluid continuum. In fact, for many purposes essential to grasp- 
ing the underlying details, the ISM is believed to possess a self- 
simil ar hierarchical structure over severa l orders of magnitude i n 
scale dLarsonlll98ll ; lFalgarone et al.lll992l ; lHeithausen et alj|l998l) . 
Direct H I absorption observations and interstellar scintillation 
measu rements suggest that the structure extends down to a scale of 
10 au dCrovisier et al]l985l ; lLan~ger et al]l995l;|Faison et al.l 19981) 
and possibly even to sub-au scales dHill etalj|2005ir Numerous 
theories have attempted to explain the origin, evolution and mass 
distribution o f these clouds and it has been established, from both 
observations ( Elmegreen &Falgaronel 19961) and numerical simula- 
tions l|Buiierteta^ 

2000), that the interstellar medium has a clumpy hierarchical self- 
similar structure with a fractal dimension in three-dimensional 
space. The main reason for this is still not properly understood, but 
it can be the consequence of an underlying fractal geometry that 
may arise due to turbulent processes in the medium. 

A theoretical study of these astrophysical systems — either a 
fluid continuum or a fractal structure — will necessitate the appli- 
cation of the mathematical principles of nonlinear dynamics. This 
is the principal objective of this work. The physical processes in 
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a fractal mediu m have been ana lysed by fractional integration and 
differentiation dZaslavsky||2002l . and references therein). To do so, 
the fractal medium has had to be replaced by a continuous medium 
and the integrals on the network of the fractal medium has had to 
be approximated by fractional integrals dRen et alj2 003). The inter- 
pretation of fractional integr ation is connected with fractional mass 
dimension (Mandelbrot 1983). Fractional integrals can be consid- 
ered as int egrals over frac tional dimension space within a numer- 
ical factor (Tar asovll2004l) . This numerical factor has been chosen 
suitably to get the right dimension of any physical parameter and to 
derive the standard expression in the continuum limit. The isotropic 
and homogeneous nature of dimensionality has also been incorpo- 
rated properly. All of these will give a se lf-consiste nt description of 
the hydrodynamics in a fractal medium ( iRovll 20071) . 

Once the hydrodynamic equivalence has been established, it 
has then been a fairly easy exercise to model the steady fr actal flow 
like a simple first-order aut onomous dynamical system (Strogatz 
ll994l:IJordan & Smitrj[l99l) . This has made it possible to gain an 
understanding of the critical aspects of the stationary phase por- 
trait of the fractal flow, especially the behaviour of the transonic 
solution. The critical point in the phase portrait has been shown 
to be a saddle point, and the transonic solution that has to pass 
through this point has been shown to be infinitely sensitive to the 
choice of a boundary condition. While this bodes ill for the fea- 
sibility of transonicity itself within the stationary framework, all 
steady global solutions (transonic or otherwise) have been found to 
be stable under a time-dependent linearised perturbation. An inter- 
esting fact that has come to light is that the necessary mathematical 
conditions, which include an equation of motion for the dynamic 
perturbation and its relevant boundary conditions, to argue for the 
stability of the steady fractal flows, have been found to be entirely 
identical to what has been reported earlier r egarding the stability 
of continuous spherica lly symmetric inflows jPetterson et al.l 19801 : 
iTheuns &DaviJl992h . This similarity is fortunate and armed with 
this knowledge, it can be safely claimed that fractal flows are just 
as stable as continuous flows under the effects of small linearised 
perturbations. 

Having said this, one will still have to confront the fact that the 
time-dependent perturbative analysis has done nothing to indicate 
the primacy of the transonic solution, and that the steady transonic 
inflow solution would not be possible without an infinite precision 
in prescribing a proper boundary condition. This obstacle has, how- 
ever, been overcome by taking into consideration explicit dynamics 
in the flow system, and then evolving a physical flow through time, 
after having started with appropriate initial conditions. Transonicity 
becomes evident very soon, and it has been argued with a simpli- 
fied analytical model in the "pressure-free" limit, that the guiding 
physical cr iterion to se lect the transonic solution is the one for- 
warded bv lBondil l l 19521) . i.e. the transonic solution will be chosen 
because it corresponds to a minimum energy configuration. While 
the "pressure-free" limit does not involve the fractal properties di- 
rectly, it has been demonstrated through a numerical integration of 
the dynamic flow equations of the fractal medium that transonic- 
ity is very much the favoured mode of infall in this case too. And 
the most salient result to have emerged from this numerical study 
has been that transonic features becomes more pronounced, as the 
medium is more like a fractal. 

It has already been mentioned that the perturbative treatment 
has been shown to give no clear-cut evidence to favour transonic- 
ity. Support for transonicity, however, has come indirectly from the 
perturbative angle too. The equation of motion for the propagation 
of the linearised perturbation has been shown to have subtle simi- 



lariti es with an eq uation implying the metric of an acoustic black 
hole dVisserl 19981) . This hints at the fact that matter might cross the 
sonic horizon at the greatest possible rate, i.e. transonically, just as 
matter has to cross the event horizon of a black hole maximally. 



2 THE EQUATIONS OF THE FLOW AND ITS FIXED 
POINTS 

Considering the existence of a medium that has a fractal structure 
of mass dimensionality D — 3d (with d < 1) embedded in a 3- 
dimensiona l space, th e mass enclosed in a sphere of radius r can be 
written as ( Roy 2 0071) 
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with D referring to the dimension, p to the constant density of the 
medium, and l c to a characteristic inner length of the medium that 
can take an arbitrary value in the limit d — > 1. This is the scale be- 
low which the medium will be continuous. The fractional integrals 
are computed as integrals over fractional dimension space within a 
numerical factor. The fractional infinitesimal length for a me dium 
with isotropic mass dimension will, therefore, be given by dRovl 
120071) 
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with the constant having been chosen to derive the standard ex- 
pression in the limit d — ► 1. It is to be noted that the infinitesimal 
area and volume elements in this "fractional continuous" medium 
of mass dimension D — 3d will be different, and hence the mass 
enclo sed in a sphere of radius r for constant density p will be dRovl 
l2007h 
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In this medium the inviscid Euler equation, de scribing t he dy- 
namics of the velocity field, v, can be expressed as dRovll2007l) 
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where <j>(r) is the gravitational potential of the central accretor 
that drives the flow (with the prime denoting the spatial deriva- 
tive of the potential). This is a local conservation law and, as it 
is to be expected, this has exactly the same form as that of the 
equation for the continuous medium. In the case of stellar accre- 
tion, the flow is driven by the Newtonian potential, <j> — —GM/r. 
On the other hand, frequently in studies of black hole accretion, 
it becomes convenient to dispense completely with the rigour of 
general relativity, and instead make use of an "effective" pseudo- 
Newtonian potential that will imitate general relativistic effects 
in the Newtonian construct of space and time dPaczvnski & Wiital 
ll98fj;lNowak & Wagoned! 1 99 lMArtemova et al.lll996h . The choice 
of a particular form of the potential will not affect the general argu- 
ments overmuch. 

The pressure, p, is related to the local density, p, through a 
polytropic equation of state p = A'p 7 , in which A' is a constant, 
and 7 is the polytropic exponent, whose admissible range is given 
by 1 < 7 < 5/3. This range is restricted by the iso thermal and the 
adiabatic limits, respectively (Chandrasekhar 1939^. The evolution 
of p is described by the equation of continuity ( Roy||2007h . 
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Figure 1. Integral solutions of the stationary fractal flow, driven by the New- 
tonian potential, for D = 2.55 and 7 = 1.4. The continuous curves, A 
and W represent "accretion" and "wind", respectively. The fixed point is at 
r = r c and |M| = 1. This point is a saddle point. 
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in which a — 3d — 1. 

The flow system will, therefore, be specified by equations 13} 
and l[5}, along with the polytropic equation of state. Of particular 
interest are the steady state solutions for the case where the fractal 
medium is at rest at a large distance from the accretor. Since tran- 
sonic flows are of primary concern here, it is required that the static 
flow should evolve from v — > as r — > 00 (the outer boundary 
condition) to v > a(r) for small r, where a is the speed of sound, 
given by a 2 = dp/ dp = ^Kp 1 ' 1 . The stationary state implies 
dv/dt = dp/dt — 0. Consequently, equations (O and {5} will be 
reduced to their steady state forms as 
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respectively. It is easily seen that equations ((6} and ((7} remain in- 
variant under the transformation v — > —v, i.e. the mathematical 
problem for inflow s (v < 0) and ou tflows (v > 0) is identical in 
the stationary state 1 Choud hurill 1 9991) . 

It is in principle possible to eliminate either v or p and solve 
for the other variable as a function of r. However, adopting a 
slightly different approach, it is possible to recast equations l[6) 
and l[7jl in a combined form as 
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whose integral solutions have been shown in Fig. [7] with <f>(r) = 
—GM/r, and with the vertical axis scaled by the absolute value 
of the Mach number, M = v/a. The two continuous curves la- 
belled "A" and "W" refer to accretion and wind, respectively. The 
critical points in the flow will be derived from the standard require- 
ment that the flow solutions will have a finite gradient when they 
will cross the sonic horizon (where the bulk flow velocity exactly 
matches the speed of acoustic propagation), which will mean that 
both the numerator and the denominator of equation l[8) will have t o 
vanish simultaneously and non-trivially ( Chakrabart3 ll990l 1 1996b . 
This can only happen when 
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which gives the critical point (or the sonic point in this particular 
case) conditions, with the subscript "c" labelling the critical point 
values. 

It is not a difficult exercise to integrate equation ® and then 
transform the variable p in it to a with the help of the equation 
of state. This, with the critical point conditions as given by equa- 
tion J9J, will give a relation for fixing the critical point coordinates 
in terms of the flow parameters, £ (which is actually Bernoulli's 
constant), a and 7 as 
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c4>'(r c ) 
2a 
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The form of 4>{r) will obviously determine the number of the criti- 
cal points, and for the Newtonian potential only one root of r c will 
be obtained from equation dlOl ). This root, for £ fixed by the outer 
boundary condition of the transonic inflow solution, will be given 
as 
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with aoo being the speed of s ound at the outer boundary of the 
flow ( Chakrabart3 [T990L Il996h . where the influence of gravity is 
negligibly weak. 

It should be worth mentioning here that although with the 
choice of a pseudo-Newtonian potential multiple roots for r c would 
result, practically speaking only one of these roots would be a phys- 
ically meaningful critical point, through which an integral solu- 
tion may pass. For spherically symmetric flows in t he general rel- 
ativistic framework, this issue has been discussed by Man dal et al.l 
d2007h . 



3 THE FLOW AS AN AUTONOMOUS DYNAMICAL 
SYSTEM 

So far the flow variables have been ascertained only at the critical 
points. Since the flow equations are in general nonlinear differen- 
tial equations, short of carrying out a numerical integration, there 
is no completely rigorous analytical prescription for solving these 
differential equations to determine the global nature of the flow 
variables. Nevertheless, some analytical headway could be made 
after all by taking advantage of the fact that equation {8j, which 



gives a complete description of the r 



phase portrait of the 



flow, is an autonomous first-order differential equation, and as such, 
could easily be recast into the mathematical form x — X(x, y) and 
y — Y(x,y), whic h is that of the very familiar coupled fir st-order 
dynamical system (Strogatzl ll994l ; Uordan & Smith! fl999). Quite 
frequently for any nonlinear physical system, a linearised analytical 
study of the properties of the fixed points of a first-order dynamical 
system, affords a robust platform for carrying out an investigation 
to understand the global behaviour of integral solutions in the phase 
portrait. 

And so it is that to investigate the nature of the critical point, 
equation ^ will have to be decomposed in terms of a mathematical 
parameter, r, to read as 
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in both of which the parameter r does not make an explicit ap- 
pearance in the right hand side, something of an especial advan- 
tage that derives from working with autonomous system s. This kind 
of pa rametrization is quite common in fluid dynamics jBohr et ail 
1 1993b . and in accretion studies especially, this approach has been 
made before (|Rav & Bhattacharieel |2002|; lAfshordi & Paczvnski 
2003 ; IChaudhurv et alj|2006t iMandal et alj|2007l : iGoswami et al 
2007). Some earlier works in accretion had also made use of the 



general mathematical aspects of this approach jMatsumoto et al.l 
1 1984 iMuchotrzeb-Czernvl 1 19861 : lAbramowicz & Katol Il989h . A 
further point that has to be noted is that the function a 2 in the right 
hand side of equation d!2b can be expressed entirely in terms of 
v 2 and r, with the help of the equation of state and equation 10. 
This will exactly satisfy the criterion of a first-order autonomous 
dynamical system. 

The next task would be to make a linearised approximation 
about the fixed point coordinates and extract a linear dynami- 
cal system out of equations d!2b . This will give a direct way to 
establish the nature of the critical points (or fixed points). Ex- 
panding about the fixed point values, a perturbation of the kind 
v 2 = v 2 + Sv 2 = v 2 (l + ei) and r = r c + Sr — r c (l + £2) 
can be applied. Using the continuity equation and the equation 
of state, this perturbation scheme, when linearised, will also give 
5a/a 2 — — (7 — 1) (ei + 2ae^) /2. Applying this perturbative 
expansion on equation d 1 2b . and linearising in t\ and £2 will give, 
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Using solutions of the type £1 ~ exp(f2r) and £2 ~ exp(S7r) 
in equations d!3b . the eigenvalues of the stability matrix associated 
with the critical points will be derived as 
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with a c itself being a function of r c , as given by equation (9). 

Once the position of a critical point, r c , has become known 
from equation dlOb , it is then quite easy to determine the nature of 
that critical point by using r c in equation d 14b . Since r c is a func- 
tion of £ and 7, it effectively implies that Q 2 can, in principle, 
be regarded as a function of the flow parameters. From the form 
of fi 2 in equation d!4b . a generic conclusion that can be immedi- 
ately drawn is that any critical point, as it may be expected for a 
conservative system, will be either a saddle point (for Q, 2 > 0) 
or a centre-type point (for Q 2 < 0). For the particular case of the 
Newtonian potential, <j> = —GM/r, the eigenvalues of the stability 
matrix will be given by 



Q 2 N = aat [(2a + 1) - 7 (2a - 1) 



(15) 



For the values of 7 and a lying in the range of physical interest, it 
can always be shown that Q. 2 ^ > 0. Hence in this case the critical 
point is a saddle point, and the curves which have been labelled "ac- 
cretion" and "wind" in Fig.[T]are in fact separatrices of a dynamical 
system, rather than physical solutions. 

The understand the full import of this line of reasoning, what 
has to be borne in mind is that saddle points are inherently unsta- 



ble, and to make a solution pass through such a point, after start- 
ing from an outer boundary condition, will e ntail an infinitely pre- 
cise fi ne-tuning of that boundary condition ( Ray & Bhattachariee 
120021) . This can be demonstrated through simple arguments. Going 
back to equations J 1 3b . the coupled set of linear differential equa- 
tions in ei and £2 can be set down as 



dei _ dei/dr _ Qi£i + Q 2 e 2 
de 2 de 2 /dr Q 3 £i + Q4£2 ' 



(16) 



in which the constant coefficients Qi, Q2, Q3 and Q4 are to be 
determined simply by an inspection of equations J 1 3b . It is also 
to be easily seen that Qi = — Q4. This makes the integration of 
equation d!6b a straightforward exercise, and it yields 



Q2S2 + 2Q 1 e 1 e 2 - Q 3 e?+C 







(17) 



with C being an integration constant. Generally speaking equa- 
tion d!7b is the equation of a conic section in the £2 — £1 plane. If 
the origin of this plane were to be considered to have been shifted 
to the saddle point, then the condition for solutions passing through 
the origin, i.e. £1 = £2 = 0, would be C = 0, which would re- 
duce equation d!7b to a pair of straight lines intersecting each other 
through the origin itself. All other solutions in the vicinity of the 
origin will, therefore, be hyperbolic in nature, a fact that is given by 
the condition (Q1 + Q2Q3) > 0. For the case of <j> = -GM/r,this 
contention can be verified completely analytically, and this shows 
that even a very minute deviation from a precise boundary condi- 
tion for transonicity (i.e. a boundary condition that will generate 
solutions to pass only through the origin, £1 = £2 = 0) will take 
the stationary solution far away from a transonic state. This ex- 
treme sensitivity of transonic solutions to boundary conditions is 
entirely in keeping with the nature of saddle points. It may be imag- 
ined that in a proper astrophysical system such precise fulfillment 
of a boundary condition will make the transonic solution well-nigh 
physically non-realisable. Indeed, this difficulty, for any kind of ac- 
creting system, is readily appreciated by anyone trying to carry out 
a numerical integration of equation l[8j to generate the transonic 
solutions, which can only be obtained when the numerics is first 
biased in favour of transonicity by using the saddle point condition 
itself as the boundary condition for numerical integration. 

Apart from this, there is also a mathematical aspect of the 
physical non-realisability of transonic solutions. Using the condi- 
tion C = will make it easy to express e\ in terms of £2 and vice 
versa. Going back to the set of linear equations given by equa- 
tions d 1 3b and choosing the second one of the two equations (the 
choice of the first would also have led to the same result), one gets 



d£2 

d7 



QI + Q2Q3) £2, 



(18) 



which can be integrated for both the roots from an arbitrary initial 
value of £2 = £2 lying anywhere on the transonic solution, to a 
point £2 = A, with A being very close to the critical point given 
by £2 = 0. Using the equivalence that Q 2 — Q 2 + Q2Q3, it can be 
shown that 



I f A 



d£2 . 1 . A 

£2 1 1 £? 



(19) 



from which it is easy to see that |r| — > 00 for A — > 0. This 
implies that the critical point may be reached along either of the 
separatrices, only after |r| has become infinitely large. This diver- 
gence of the parameter r indicates that in the stationary regime, 
solutions passing through the saddle point are not actual solutions, 
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but separatrices of v arious classes of solutions JStrogatz|[l99i 
I Jordan & Smithl 19991) . This fact, coupled with the sensitivity of the 
stationary transonic solutions to the choice of an outer boundary 
condition, makes their feasibility a seriously questionable matter. 



4 A TIME-DEPENDENT PERTURBATIVE APPROACH 

It has been demonstrated in the foregoing section that the steady 
transonic accretion solution is unstable under infinitesimal devia- 
tions from the precise outer boundary condition needed to generate 
the solution. In the astrophysical context, such precision is quite im- 
possible, and, therefore, the very feasibility of transonicity becomes 
a matter of grave doubt. This difficulty can, however, be dispelled 
if one is mindful of the fact that the real astrophysical flow is not 
static but dynamic in character, i.e. it will have an explicit depen- 
dence on time. This, of course, will mean that the time-dependent 
terms involving both the velocity and the density fields in equa- 
tions © and (|5) will have to be retained. 

While full-fledged time-dependence of the flow variables will 
undoubtedly reveal many new interesting mathematical facets (all 
of them involving the mathematics of partial differential equations) 
of the physical problem, it would still be worthwhile to go back to 
studying the properties of the background stationary flow under the 
influence of a linearised perturbative effect. As a preliminary ex- 
ercise in accounting for explicit time-dependence, this will at least 
shed some light on the global stability of the flow solutions. 

To that end, it will first be necessary to define, closely f ollow- 
ing a perturbative procedur e prescribed by IPetterson et alj dl980h 
and lTheuns & David Jl992h . a new physical variable / = pvr a . It 
is quite obvious from the form of equation l[5) that the stationary 
value of / will be a global constant, fo, which can be closely iden- 
tified with the matter flux rate. A perturbation prescription of the 
form v(r, t) = vo(r) + v'(r, t) and p(r, t) = po(r) + p'(r, t), will 
give, on linearising in the primed quantities, 



\po v 



(20) 



with the subscript "0" denoting stationary values in all cases. From 
equation ((5), it then becomes possible to set down the density fluc- 
tuations, p', in terms of /' as 



dp' vppo 
dt fo 



df_ 

dr 



0. 



(21) 



Combining equations J20b and d21t will then render the velocity 
fluctuations as 

dv' 
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' df df 

;» ~ f V dt +Vo dr 
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which, upon a further partial differentiation with respect to time 
will give 

dt 2 



d 
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v dr ). 
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From equation © the linearised fluctuating part could be ex- 
tracted as 



dv' 8 ( , 2 p' \ 
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with ao being the speed of sound in the steady state. Differentiating 
equation d24t partially with respect to t, and making use of equa- 
tions ( 121b . ( 122b and ((23} to substitute for all the first and second- 
order derivatives of v' and p', will deliver the result 
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A little readjustment of terms in equation ( 1251 ) will finally give an 
equation of motion for the perturbation as 

d 2 f n d ( df'\ Id 
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which is an expression that is exactly the same as what can be 
derived upon perturbing the stationary sol utions of spherically 
symmetric inflows in a continuous medium JPetterson et alj[l98u 
Theuns & David 1992). Another aspect of equation ([26) is that its 
form has no explicit dependence on the potential — Newtonian or 
pseudo-Newtonian — that is driving the flow. This is entirely to 
be expected, because the potential, being independent of time, will 
only lend its direct presence to the stationary background flow. Ar- 
guments regarding stability will, therefore, be more dependent on 
the boundary conditions of the steady flow. As the form of the equa- 
tion of motion for the linearised perturbation remains unchanged 
even for a flow in a fractal medium, and as the physical boundary 
conditions are al so not altered in this case, the general conclusions 
reached by both lPetterson et alj dl980h and lTheuns & David ( 1 1992b 
regarding flows in a continuous medium, will carry over here, and 
it can be safely claimed that under all reasonable boundary condi- 
tions, both the transonic and subsonic solutions will be stable. 

While this does nothing to cause any immediate worry, it also 
does not reveal anything in particular either about the physical fea- 
sibility of any solution from a perturbative point of view. This is in 
keeping with the conve ntional wisdom about spherically symmet- 
ric flows lGarlick] [l979h that the natural preference of the system 
for any particular solution — especially the transonic solution — 
cannot be justified by a linear stabilit y analysis, b ut by the more 
fundamental arguments forwarded bv lBondil ( 119521) . 

For all that, some positive hint about the primary status of the 
transonic solution can actually be derived if the whole issue of a 
linear stability analysis is viewed from a different perspective. It 
is known that there is a close one-to-one correspondence between 
certain features of black hole physics and the physics of supersonic 
acoustic flows dVissedll998h . In this very context, a compact ren- 
dering of equation ((25) can be obtained as 



9m (t^duf) = 0, 



(27) 



in which the Greek indices are made to run from to 1, with the 
identification that stands for t, and 1 stands for r. An inspection 
of the terms in the left hand side of equation d25t will then allow 
for constructing the symmetric matrix 



Vo 

To 



Vo 



a 2 , 



(28) 



Now in Lorentzian geometry the d' Alembertian for a scalar field in 
curved space is given in terms of the metric g M „ by ( Visser 1998) 

1 



Aip 



--d u 



(29) 



with being the inverse of the matrix implied by g^. Compar- 
ing equation ((27) with equation ((29), it would be tempting to look 
for an exact equivalence between f" and -^/—gg' ly . This, how- 
ever, cannot be done in a general sense. What can be appreciated, 
nevertheless, is that equation 027 1 gives a relation for /' which is of 
the type given by equation {29). The metrical part of equation ( 127b . 
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as given by equation l !28t . may then be extracted, and its inverse 
will incorporate the notion of a sonic horizon of an acoustic black 
hole when Vq = 0%. This point of view does not make for a per- 
fect acoustic analogue model, but it has some similar features to the 
metric of a wave equation for a scalar field in curved space-time, 
obtained through a somewhat different approach, in which the ve- 
locity of an an irrotational, inviscid and barotropic fluid flow is first 
represented as the gradient of a scalar function if), i.e. v = — V-0, 
and t hen a perturbation is imposed on this scalar function jVissed 
1 1998b . 

The foregoing discussion indicates that the physics of super- 
sonic acoustic flows closely corresponds to many features of black 
hole physics. This closeness of form is very intriguing. For a black 
hole, infalling matter crosses the event horizon maximally, i.e. at 
the greatest possible speed. By analogy the same thing may be said 
of matter crossing the sonic horizon of a spherically symmetric 
fluid flow, falling on to a point sink. That this fact can be appre- 
ciated for the spherically symmetric accretion problem, through a 
perturbative result as given by equation l |25t , is quite remarkable. 
This is because it is universally recognised that that no insight into 
the special status of any inflow solution may possibly be derived 
solely through a perturbative technique dGarlickl 1 1 979T) . It is the 
transonic s olution that c rosses the sonic horizon at the greatest pos- 
sible rate (Bondi 1952), and the similarity of form between equa- 
tions l |25t and l !29t may very well be indicative of the primacy of 
the transonic solution. If such an insight were truly to be had with 
the help of the perturbation equation, then the perturbative linear 
stability analysis might not have been carried out in vain after all. 



5 DYNAMIC EVOLUTION TOWARDS TRANSONICITY 

Much more direct and robust evidence in favour of transonicity 
could be obtained if the accreting system were to be made to evolve 
through time, as opposed to making it suffer small linearised pertur- 
bations in time. Having said this, it must also be stressed that equa- 
tions © and l[5}, which govern the temporal evolution of the flow, 
are not amenable to a ready mathematical analysis; indeed, in the 
matter of incorporating both the dynamic and the pressure effects 
in the equations, short of a direct numerical treatment, the math- 
emati cal problem, is very aptly described as "insuperable" dBondil 
1952). Therefore, to have a preliminary appreciation of the govern- 
ing mechanism that underlies any possible selection of a transonic 
flow, it should be necessary to adopt some simplifications. This will 
pave the way for a more complete physical understanding of the 
evolutionary properties of the flow. 

The evolutionary dynamics is, therefore, to be studied first in 
the regime of what is understood to be the "pressureless" motion 
of a fluid in a gravitational field Jshulll99lh . as opposed to drop- 
ping t he dynamic effects to study a much simplified stationary pic- 
ture dBondil [I952I) . Simplification of the mathematical equations, 
however, is not the only justification for such a prescription. A 
greater justification li es in the fact t hat the result delivered is in 
conformity wi th, whatlGarlickl d 19791) calls "the more fundamental 
arguments' 'of lBondil ( fl95lk that it is the criterion of minimum to- 
tal energy associated with a solution, that will accord it a principal 
status over all the others. 

An immediate consequence of adopting dynamic equations is 
that the invariance of the stationary solutions under the transforma- 
tion v — > — v, is lost. As a result, one will have to separately con- 
sider either the inflows (v < 0) or the outflows (v > 0), a choice 
that has to be imposed upon the system at t = 0. Euler's equation, 



tailored according to the simplified requirements of a "pressure- 
less" field, is rendered as 



dv ^dv GM _ 
dt dr r 2 



(30) 



which can be solved by the method of characteristics jDebnathl 
1 19971) . The characteristic solutions are obtained from 



dt 
1 



dr 

v 



dv 



(3D 



(32) 



-GM/r 2 ' 

First solving the dv / dr equation will give 

tP_ _ GM_ _ <?_ 
Y r~ ~ ~2 
in which c is an integration constant that derives from the spatial 
part of the characteristic equation. This result is to be used to solve 
the dr /dt equation from equation <3U . which will finally lead to 

2vr , , /v \ 2 2ct 
In r — m — hi = c 

cr s \c J r s 



(33) 



with c being another integration constant, and r s being a length 
scale in the system, defined as r s = 2GM/(? . 

A general solution of equation d3 1 b is given by the condition, 
/(c) = c 2 /2, with / being an arbitrary function, whose form is to 
be determined from the initial condition. The general solution can, 
therefore, be set down as 

2 



GM 
r 



f 



2vr 
cr s 



In r 



m( v - 

> c 



2ct 

r s 



(34) 



to determine whose particular form, the initial condition that will 
have to be used is v = at t = for all r. This will lead to 

-.2 



GM 
r 



GM (v , 
+ 1 

r V c 



-2 2vr 2cf\ 

exp (35) 



from which it is easy to see that for t 
is the stationary solution, 



GM 



0. 



oo, what is approached 



(36) 



Corresponding to the given initial condition, this is evidently the 
stationary solution associated with the lowest possible total energy, 
and the temporal evolution selects this solution from among all the 
others. The whole picture could be conceived of as one in which 
a system with a uniform velocity distribution v = everywhere, 
suddenly has a gravity mechanism switched on in its midst at t = 0. 
This induces a potential —GM/r at all points in space. The system 
then starts evolving to restore itself to another stationary state, with 
the velocity increasing according to equation l |35l l, so that for t — > 
oo, the total energy at all points, E = (v 2 /2) - (GM/r) = 0, 
remains the same as at t = 0. 

This contention has been borne out by a numerical integration 
of equation J30t by the finite differencing technique. The mass of 
the accretor has been chosen to be A/0, while its radius is r®. The 
evolution through time has been followed at a fixed length scale of 
Sir©. The result of the numerical evolution of the velocity field, 
—v (for inflows v is actually negative), through time, t, has been 
plotted in Fig. [2] The limiting value of the velocity, as the evolu- 
tion progresses towards the long-time limit, is evidently \j2GM/r 
(with M = Mq and r = 51r ), as equation l |36| > would suggest. 
This is what the plot in Fig.|2]shows, as —v approaches its terminal 
value for t — > 00. The slope of this logarithmic plot also indicates 
that in the early stages of the evolution there is a linear growth of 
the velocity field through time, but on later times, conspicuous de- 
viation from linearity sets in. 
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Figure 2. For the temporal evolution of velocity, the slope of this logarith- 
mic plot shows that in the early stages of the evolution, — v varies linearly 
with t. Deviation from linear growth sets in later. The horizontal line on 
top shows the terminal value of the velocity field, whose evolution is being 
followed at a fixed length scale. 
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Figure 4. Velocity field, scaled as the Mach number, M, after 4000 seconds, 
for various values of d. Moving from the bottom to the top, successive so- 
lutions have been shown for d = 1.0, 0.85 and 0.7, respectively. The radial 
distance has been normalised with respect to r©, and plotted logarithmi- 
cally over a length scale of r© to 4r© . 
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Figure 3. The log-log plot of the Mach number versus the radial distance, 
with the latter being scaled by the sonic radius, r c , after 4000 seconds of 
evolution. Going from left to right each curve corresponds to d = 0.7, 
0.85 and 1.0, respectively. The slope of the curves indicates a power law 
behaviour for the evolution on intermediate length scales. All the curves 
approach a saturation velocity scale closer to the accretor. 

The argument presented above, with the effects of pressure 
taken into account, can now be extended to understand the dynamic 
selection of the transonic solution. The inclusion of the pressure 
term in the dynamic equation, fixes the total energy of the system 
accordingly at t = 0. A physically realistic initial condition should 
be that v = at t = 0, for all r, while p has some uniform value. 
The temporal evolution of the accreting system would then non- 
perturbatively select the transonic trajectory, as it is this solution 
with which is associated the least possible energy configu ration. 
This a rgument is in conformity with the assertion made by iBondJ 
that it is the criterion of minimum total energy that should 
make a particular solution (the transonic solution in this case) pre- 
ferred to all the others. However, this selection mechanism is effec- 
tive only through the temporal evolution of the flow. 

To test this contention a numerical study has been carried 
out once again using finite differencing, but this time using both 
the dynamic equations for the velocity and the density fields, as 



given by equations I0 and l[5)- The accretor has been chosen to 
have a mass, M©, and radius, r©. The "ambient" conditions are 
aoo = 10 km s _1 and = 10 -21 kg m~ 3 . The polytropic 
exponent, 7, has been set as n = (7 — = 1.61. For these 

values of the physical constants, transonicity becomes apparent 
even at the very early stages of the evolution. This is shown in 
Fig. [3] in which the velocity field (scaled as the Mach number) 
has been plotted after it has evolved for 4000 seconds. The hori- 
zontal distance has been scaled by the sonic radius (which, going 
by equation Jilt , effectively makes this scaling dependent on the 
fractal dimension), and it shows a saturation behaviour for the in- 
flow velocity at smaller length scales. This saturation scale for the 
velocity field is roughly the same for any value of d, although at 
larger length scales, a curve placed higher in the plot, corresponds 
to a higher value of d. Besides this, all curves, as it is apparent 
from the slope of each of them, display the same kind of power- 
law behaviour on larger length scales. All this is somewhat remi- 
niscent of the growth processes exhibited by the ballistic deposi- 
tion model, which is applied t o generate a nonequilibrium inter- 
face Igarabasi & Stanlevlll995l) . 

The comparative properties of the velocity field, at smaller 
length scales, for various values of d, have been exhibited in Fig. [4] 
Here the radial distance has been scaled in terms of the radius of 
the accretor, which in this case is r©. What is interesting to note 
in this plot is that for all other conditions remaining the same, on 
length scales close to the accretor, solutions corresponding to lower 
values of the fractal dimension, d, grow faster in time than solu- 
tions with higher values of d. It is possible to argue that this is 
exactly how it should be. The pressure of the infalling gas, in so 
far as it is connected to the density through a polytropic equation 
of state, builds up resistance against gravity, because of the growth 
of the density field on small length scales. Transonicity can only be 
achieved when gravity wins over pressure on length scales smaller 
than the sonic radius. This will be all the more true near the accre- 
tor, where the velocity field will evolve under free-fall conditions, 
and, therefore, the more dilute the gas, the more efficient will be 
the drive towards the transonic state. Now a fractal medium may 
be viewed equivalently as a continuum with an effective lesser den- 
sity. In this situation a system with a lesser value of d will be more 
prone to losing against gravity than a system with higher value of 
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d, and so the race towards transonicity will be more successful as 
d decreases. It is exactly this state of affairs that Fig.[4]graphically 
represents. 



6 CONCLUDING REMARKS 

An earlier work reported bv lRovl d2007h was carried out under the 
implicit assumption that the accretion process would take place 
transonically. The present treatment b ears out this assumption self- 
consistently. It was also discussed by IrovI J2007h that the rate of 
accret ion in a fractal medium can vary significantly from the lBondj 
rate for more massive accretors. This is very much in con- 
formity with the conclusions derived in this work, through the 
numerical evolution of transonicity, and it would be judicious to 
account for this fact, while studying the accretion of a fractal 
medium on to a black hole . Bla ck hole accretion is necessarily 
transonic dChakrabartilll990lll996h . but even for accretion from a 
molecular cloud on to a star, in the absence of any inner bound- 
ary condit ion being imposed on the flow, the flow is expected to be 
transonic jPetterson et al.|[l98(il) . Therefore, whatever be the nature 
of the accretor, both transonicity and the quantitative modifications 
arising due to the fractal nature of the accreting medium, will be 
very much relevant for studies in spherically symmetric accretion. 

It has also been discussed here that transonic properties man- 
ifest themselves more noticeably with an increase in the fractal 
properties of the flow. The dynamic evolution has shown that the 
growth rate of the velocity field (as scaled against the speed of 
acoustic propagation) becomes significantly higher in this case. 
This, of course, will have a direct bearing on the mass accretion 
rate, and it is worth conjecturing that there might be some observa- 
tional evidence for this kind of behaviour. 
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